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We present an estimate of the bolometric X-ray luminosity - velocity dispersion 
(L x — a v ) relation measured from a new, large and homogeneous sample of 171 low 
redshift, X-ray selected galaxy clusters. The linear fitting of log(-L x ) — log(cr„) gives 
L x = io 32 - 72±0 08 '4.i±o.3 erg s -i j 1 -2 i Furthermore, a study of 54 clusters, for which 
the X-ray temperature of the intracluster medium T is available, allows us to explore 
two other scaling relations, L x — T and a v — T. From this sample we obtain L x oc 
ji3.i±o.2 anc j a ^ x jii.oo±o.i6^ wmcn are fully consistent with the above result for the 

L x -a v . The slopes of L x — T and a v — T are incompatible with the values predicted 
by self-similarity [L x oc T 2 oc a^), thus suggesting the presence of non-gravitational 
energy sources heating up the intracluster medium, in addition to the gravitational 
collapse, in the early stages of cluster formation. On the other hand, the result on 
log(i x ) — log((T„) supports the self-similar model. 
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1 INTRODUCTION 

Galaxy clusters are dark matter haloes which form mainly 
through the process of gravitational collapse, and so they are 
expected to be self-similar: small clusters or galaxy groups 
must be scaled-down versions of the more massive systems. 
Self-similarity predicts a series of relationships between dif- 
ferent cluster observables such as X-ray luminosity L x , mass 
M, temperature T and galaxy velocity dispersion a v , which 
are the same in all systems (Kaiser 1991). 

The radial velocity dispersion a v of the galaxies in a 
cluster probes the depth and shape of the potential well, 
assuming that the luminous matter traces reasonably well 
the dark matter in clusters. Moreover, the intracluster gas 
emits X rays through a process of thermal bremsstrahlung 
and its bolometric luminosity L x is found to be strongly 
correlated with a v , as both gas and galaxies share a common 
potential well (Solinger & Tucker 1972). Also, since the X- 
ray emission can be well modelled by thermal emission of a 
hot, optically thin plasma, L x and the plasma temperature T 
must be correlated. Finally, because galaxies are embedded 



* Partially based on observations at the European Southern Ob- 
servatory La Silla, Chile 



in the intracluster medium, a correlation is also expected 
between the gas temperature T and the galaxies velocity 
dispersion a v . 

Self-similar models consider that the only energy source 
in the cluster comes from the gravitational collapse, predict- 
ing the following scaling relations: L x oc T 2 oc cr*. Whereas 
there seems to be a general agreement between different 
measurements that L x oc T~ 3 , the measurement of L x — a v 
has given somewhat contradictory results so far. Some au- 
thors have found that L x oc a 4 indeed, although with quite 
large measurement errors or rather small data samples, while 
others find slopes larger than 4 (see Table 0. Part of the 
differences in the results could come from different ways of 
selecting the sample, with a preference for regular clusters 
in some of these surveys. 

It has also been suggested that clusters and groups do 
not follow the same L x — a v scaling relation, the latter being 
flatter than the former (e.g. Mahdavi et al. 2000, Xue & Wu 

2000, Dell' Antonio et al. 1994). But there are other mea- 
surements contrary to that conclusion (Mahdavi & Geller 

2001, Mulchaey & Zabludoff 1998). 

For more distant clusters (z between ~ 0.15 and ~ 0.6) 
there is some evidence that the slope is also > 4 (Borgani 
et al. 1999, Girardi & Mezetti 2001), although only small 
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samples are available at the moment, and more data are 
needed to reduce the error bars. 

There is a clear need of new measurements performed on 
large samples (more than 100 clusters), selected from a ho- 
mogeneous dataset. The largest cluster (not group) samples 
available up to date, on which this kind of study has been 
performed, are literature compilations combining data from 
different authors and/or instruments (for example Xue & 
Wu 2000, Mahdavi & Geller 2001). The present work comes 
to partially fill this gap, as we have a large homogeneous 
subsample of clusters selected from REFLEX, a flux-limited 
X-ray selected cluster catalogue built under well defined se- 
lection criteria. 

The structure of this paper is as follows. In Section[5]we 
give a brief introduction to the REFLEX catalogue and the 
way in which the subsample in this paper has been selected. 
Section |21 illustrates in detail the process of the L x and a v 
measurements, which lead to the the fitting of the log (La,) — 
log(<T„) in Section 0] In Section f5.1l we propose the removal 
of those clusters in the sample which are multiple systems. 
In Section [5] we analyse possible biases due to the nature 
of the selected sample. A study on the L x — T and a v — T 
relation follows in Section [U and we finish with a summary 
in Section|7| We are assuming Ho = 50 km s" 1 Mpc" 1 in a 
flat universe with S1a = 0.7. 



2 THE DATA 

The REFLEX (ROSAT-ESO Flux-Limited X-ray) cluster 
survey (Bohringer et al 2001, Bohringer et al 2003 in prepa- 
ration, Guzzo et al. 2003 in preparation) has identified and 
measured successfully the redshift for all southern galaxy 
clusters (at galactic latitude \bn\ > 20°) down to a flux 
limit f x > 3 ■ 10~ 12 erg s" 1 cm" 2 in the ROSAT All-Sky 
Survey. This sample has already provided us with results 
concerning the measurement of the cluster X-ray luminosity 
function at low redshift (Bohringer et al 2002), the power 
spectrum (Schuecker et al. 2001a), the spatial correlation 
function (Collins et al. 2000), the measurement of f2 m and 
us (Schuecker et al. 2003a) or the measurement of the equa- 
tion of state parameter w of the dark energy (Schuecker et al. 
2003b), among others. A detailed description of the sample 
construction can be found in Bohringer et al. (2001). 

REFLEX is an X-ray selected sample and the first com- 
pilation of the REFLEX catalogue includes 452 clusters 
down to a flux limit of 3 x 10~ 12 erg s" 1 cm" 2 in the ROSAT 
energy band from 0.1 to 2.4 keV. The survey covers a to- 
tal area in the sky of 4.24 ster. Two regions corresponding 
to the Magellanic clouds and the Galactic plane have been 
excluded because of the difficulties in identifying clusters of 
galaxies in these densely crowded regions, as well as some 
strips of the sky which were not observed by ROSAT due to 
power-off of the detector when the satellite was crossing the 
Earth's radiation belts (see Figure 2 in Bohringer et al 2001 
for a survey exposure map). 

Optical follow-up has been performed through long-slit 
and multi-object spectroscopy at the 1.5m, 2.2m and 3.6m 
ESO telescopes at La Silla (Chile) during a series of observ- 
ing campaigns, which started back in 1991. We complement 
our sample of galaxy redshifts by looking into the literature 
(using the NASA/IPAC Extragalactic Database, NED) and 



adding those redshifts already published for the previously 
known clusters. 

The cluster membership of each individual galaxy has 
been decided based on the search for velocity peaks in the 
velocity distribution of the galaxies observed in the field. 
All galaxies were assigned the same weight in the computa- 
tion of the cluster's median redshift because in most cases 
not enough photometric material is available to give proper 
weights (i.e. to identify a real cD galaxy). Anyway, Bright 
Cluster Galaxies (if present) were used to decide on the 
main redshift of the system when it showed multiple veloc- 
ity peaks in the redshift histogram along the line of sight. In 
other cases the choice among multiple velocity peaks was de- 
cided taking into account the correlation between the peak 
position and the X-ray emission centre and also the strength 
of the peak. Finally, other clusters ended up being split into 
different systems when multiple X-ray emission centres as- 
sociated to them were detected. 

We have selected from REFLEX those clusters with 
more than seven measured galaxy redshifts (see Fig. 0, in- 
cluding only galaxies at a maximum projected distance from 
the cluster center of 0.5 h^ 1 Mpc, due to the lack of data at 
larger distances from the cluster centre. Outliers in the ve- 
locity space were rejected by applying the 3 sigma clipping 
method: galaxies with velocity differences with the cluster 
centre velocity larger than three times the measured veloc- 
ity dispersion were not used in the subsequent iterations 
to avoid possible contamination by foreground/background 
objects. By following this procedure we end up with 171 
clusters out of the 452 in the whole catalogue. It is a low- 
redshift sample, the median redshift being 0.076. Only 14 
clusters have z above 0.2 (see Fig. |2J. 



3 MEASUREMENT OF L x AND ay 

3.1 Bolometric X-ray luminosity L x 

Counts in the ROSAT images were measured using the 
Growth Curve Analysis (GCA) and then converted into 
fluxes, as described in Bohringer et al. (2001). Luminosi- 
ties have first been transformed into the ROSAT rest-frame 
band (0.1-2.4 keV) (see Bohringer et al. 2002 for details) 
and afterwards into bolometric luminosities. We assumed a 
MEKAL spectral model with an ICM metallicity of 0.3 Z Q 
(Anders & Grevesse 1989). We used XSPEC (Arnaud 1996) 
to obtain the relation between the luminosity in the 0.1- 
2.4 keV band and the bolometric one, as a function of the 
cluster temperature(see Bohringer et al. 2003 for a detailed 
conversion table), which was computed from the velocity 
dispersions assuming the empirical relation between a v and 
T measured by Girardi et al. (1996). 



3.2 Radial velocity dispersion a v 

The radial velocity dispersions a v have been computed by 
means of the biweight estimator of central location and scale 
(within ROSTAT by Beers, Flynn & Gebhardt 1990). This 
estimator was first suggested by Tuckey (in Andrews et al. 
1927) as a better way to study non-gaussian or contam- 
inated normal distributions than the gaussian estimators 
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Figure 2. Rodshift distribution of the 171 clusters in the selected 
sample from the REFLEX catalogue. 



Figure 1. Distribution of the number of redshifts per cluster in 
the 171 subsample of the REFLEX catalogue which we study in 
this paper. The dashed vertical line shows the limit of 7 redshifts 
per cluster imposed to select the clusters for which a a v is mea- 
sured. 



(mean and standard deviation). The biweight location es- 
timator belongs to the family of estimators known as M 
estimators of location. It works by minimizing a function 
of the deviations of each observation from the estimate of 
location and has the advantage that it does not assume an 
underlying gaussian distribution, which in our case means to 
have a perfectly relaxed cluster and might not be the case 
in most of the sample. It is less affected by points in the 
wings of the distribution and is robust for a broad range of 
non-gaussian underlying populations. 

We used the biweight estimator of central location and 
the gapper estimator of scale (based on the gaps between 
order statistics, Wainer & Thissen, 1976) when the number 
of redshifts available in the cluster was 7 ^ N z ^ 10. For 
N z > 10 the biweight estimator was chosen for both loca- 
tion and scale. Errors were obtained in all cases by jacknifing 
of the biweight. The choice of the different estimators has 
been done on the basis of the tests carried out by Beers et al. 
(1990) for samples with different number of data points. The 
usual cosmological correction and the correction for velocity 
errors (Danese, de Zotti & Tullio, 1980) were also applied. 
The whole procedure was performed iteratively until the re- 
sults converged. 

The clusters with the larger errors in the velocity dis- 
persion estimate are among the ones for which N z ^ 10. 

One source of concern is the fact that we are computing 
the a v using galaxies which are quite close to the cluster 
centre, at distances less than 0.5 h^, 1 Mpc (median distance 
is 0.04^ Mpc). Girardi et al. (1996) find that the value of 
a v is a function of the distance to the cluster centre, and it 
reaches a stable value when galaxies at a projected distance 



of ~ 1 h^ 1 Mpc are considered. At smaller distances, the <j v 
does not show a constant behaviour: sometimes it raises up, 
sometimes it goes down, and in some cases it remains stable. 
So, in a large sample like ours we expect that this effect 
will cancel out. From our current data, it is not possible to 
determine the a v at clustercentric distances larger than 0.5 
h^ 1 Mpc. 



4 THE L x (Tv RELATION 

A power law was fit to the \og(L x ) — log(a„) relation 
by means of the Orthogonal Distance Regression method 
(ODR, Boggs et al. 1992), which takes into account errors 
on both variables. We took the data in units of L^/10 45 erg 
s -1 (L45) and velocity dispersions in units of <r„/500 km s~ x 
( o"5oo ), that is, we put the origin of coordinates more or less 
at the center of the datapoint cloud to help the model find 
the right Y-axis intercept. We will use these units through- 
out the paper, but for Tabled where we quote L x in erg s~ 
and a v in km s^ 1 for comparison purposes with the works 
referenced there. The best fit to log(L4s) — log(asoo) is found 
to be (see Figure [3J 

log(L 45 ) = (-1.34±0.08)+(4.1±0.3)-log((75oo) ergs" 1 !^^) 

The Spearman rank correlation coefficient for this sam- 
ple is r 3 = 0.51 and the probability that a random distribu- 
tion had this value of r s or larger by chance is P ^ 2.6 10~ . 

The distribution of errors in the sample has been as- 
sumed to be normal in the ODR fitting. To test for this as- 
sumption, we have performed a maximum likelihood analysis 
with sigma clipping to find which is the best fitting gaus- 
sian to the actual error distribution in X. As errors in X are 
far larger than errors in Y, we can neglect the effect of the 
latter. We find that the best gaussian has central location 
x — 0.04 and standard deviation a = 1.9 (Figure 0J. The 
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Figure 3. Bolometric X-ray luminosity L x vs. radial velocity 
dispersion a v for a REFLEX sample of 171 clusters. The linear 
fit gives a log-slope of 4.1 ± 0.3. and the dashed lines are the 
1-sigma errors. Filled circles correspond to clusters at redshift 
z sC 0.05, open triangles are clusters with 0.05 < z $C 0.1 and 
open circles are clusters at z > 0.1. 



Figure 4. Distribution of residuals along the X axis in the 171 
cluster sample. We also show the corresponding maximum like- 
lihood normal distribution. The distribution of residuals is well 
approximated by a gaussian distribution, with central location at 
x = 0.04 and standard deviation a = 1.9. 



error distribution along the X axis is thus reasonably gaus- 
sian. 

We can also perform a crude estimate of the intrinsic 
dispersion of the sample, £i n t, 



Emi = yj {a 2 b 



(2) 



in error units. a b s is the best fitting gaussian a parameter 
as found above for the current error distribution, that is 
fobs = 1-9. (Jmeas is the a parameter of the error distribution 
in the ideal case in which only the measurement errors are 
responsible for the dispersion in the data, which would be 
a gaussian with a = 1. We find Ei nr = 1.62 in error units, 
that is, the intrinsic dispersion of the sample is slightly larger 
than the dispersion due to the errors in the data. To see what 
this means in terms of physical units, let us assume that a 
typical value of the error in the velocity dispersion is the 
median of the error distribution, T, typ = 120 km s _1 . The 
intrinsic dispersion is then 



(Tint 



= 1.62 ■ St^p = 195 km s" 



(3) 



Finally, to test for the stability of the fitting, we re- 
peated it by removing those points with an error in the ve- 
locity dispersion larger than 30 percent of their value. We 
find 



log(L 4 



-1.38 ± 0.09) + (4.2 ± 0.3) • log(cr 5 



(4) 



thus showing the robustness of the fitting against the pres- 
ence of large errors in a v . 

We may also worry about the small number of redshifts 
(7 at minimum) with which we are computing the a v in 



some clusters. We have repeated the fit,this time taking only 
clusters for which a v had been determined from at least 30 
individual redshifts. In this way, we are left with a sample 
of only 57 clusters. The best fit relation is 



log(L 45 ) 



-1.44 ± 0.12) + (4.2 ± 0.4) • log (cr 5 



(•») 



which is in agreement at the one sigma confidence level with 
the fit obtained for the whole sample. Figure |S] shows the 
best fit found. 

The best fitting gaussian to the error distribution in X 
is centered at x = 0.3 and has standard deviation a = 1.9. 
This gives a value for the intrinsic dispersion of E; nt = 195 
km s _1 . 



5 DISCUSSION 

Self-similarity in the purely gravitational hierarchical sce- 
nario of structure formation implies that for galaxy clusters 
L x oc T 2 oc cr%. Although several authors have found a good 
agreement between this prediction and their measurements 
for the L x — a v relation, others have found steeper values 
(see Table GJ. 

Our value of a = 4.1 ± 0.3 is in good agree- 
ment with the slope measu re d by lOjiintari^^^Melnir^k 

il99a) . iMahdavi fc Gellei 
ll200l|). As for the intercepts , 
( 199, 



|l982[)lMulchaev fc Zabludofj 
l|200 lli and lGirardi fc MezzettJ 



our r esult is compatible with Mulchac v fc Zabludofj w .. 
and IMahdavi fc Gelled fcOOll) at about the one-sigma 
confidence level. 

From Fig.|S|one can see that our sample is mostly pop- 
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Table 1. Compilation of literature values obtained for a and b in the relation log(L^) = a + b X log(<Ti,). The number of clusters that 
entered the computation is also shown, together with some information about the sample and the cluster selection method. L x is in 
units of erg s — 1 and a v is in km s" 1 . All of them, except for Quintana & Melnick 1982, make use of bolometric luminosites. Only White 
et al. 1997 remove cooling flow clusters from the sample. Multicomponent clusters are explicitely removed in Girardi & Mezzetti 2001, 
Borgani et al. 1999 and in this work(b).This work(c) is the result obtained when using a volume-limited sample. 



Reference b a No. clusters Flux-limited sample Selection method 
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31.61 ± 1.09 


38 


no 


optical & X-ray 
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Xue & Wu (20001 
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51 


yes 


X-ray 




1 

a 7 (500 km s ') 

Figure 5. Bolometric X-ray luminosity L x vs. radial velocity 
dispersion a v for a REFLEX sample of 57 clusters for which more 
than 30 galaxy redshifts entered the a v computation. The linear 
fit gives a log-slope of 4.1 ± 0.4. and the dashed lines are the 
1-sigma errors. Filled circles correspond to clusters at redshift 
z ^ 0.05, open triangles are clusters with 0.05 < z sC 0.1 and 
open circles are clusters at z > 0.1. 



ulated by clusters and groups are practically absent, if we 
consider as "groups" those systems with a v < 340 km s~ 
(following the criterium used by Mahdavi et al. 2000) . More 
specifically, only 11 out of the 171 systems are groups. Some 
authors (e.g. Helsdon & Ponman 2000, Mahdavi et al. 2000, 
Xue & Wu 2000, Dell'Antonio et al. 1994) have found that 
the scaling laws are different for clusters and groups. In par- 
ticular, they find that clusters show a steeper slope in the 
L x — <Jv relation. Unfortunately, e cannot derive a meaning- 
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Figure 6. Velocity dispersion distribution in the cluster sample 



ful L x — a v relation for groups from our sample, due to the 
small number of groups observed. 

The complexity involved in the measurements of a v and 
L x (because of the presence of substructures, cooling flows, 
and other physical phenomena) , the largely unknown intrin- 
sic dispersion of their relationship, and the use of different 
fitting methods could explain the dif ferent results reac hed 
by the different works. For example, IXue fc Wul i200Jl il- 
lustrate the different results obtained by using ODR regres- 
sion and Ordinary Least Squares (OLS) regression in the 
analysis. 

5.1 A sample free from clusters with 
substructures. 

Substructures in clusters constitute a problem in the radial 
velocity dispersion estimates, as the different components 
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show different a v and an estimate of a single a v for the clus- 
ter may not be meaningful. We cannot make a reliable study 
of the substructures in our clusters as in many cases we do 
not have enough redshifts to perform it. The sensitivity for 
substructure detection depends strongly on the number of 
galaxies with measured redshift available for each cluster. In 
our sample, though, we have large differences in the number 
of galaxies per cluster and trying to detect substructure is 
not feasible. 

In a paper by Schuecker et al (2001b) they identify the 
clusters in the REFLEX sample which show multiple com- 
ponents in the X-ray emission. When we take those clusters 
out from our sample, we end up with 123 clusters. In this 
case we have 



log(L 45 ) = (-1-28 ± 0.10) + (4.2 ± 0.4) • log^soo) 



(6) 



which is again in very good agreement with the one 
obtained for the whole sample. 

This result shows that our sample is not largely affected 
by the problems of substructure, probably because we are 
selecting galaxies very close to the X-ray peak emission, that 
is, to the central potential well, where one would expect 
that the strength of the gravitational potential will erase 
any substructures which might be present. Even in the case 
of redshifts taken from the literature (see Section l-"i.2> we 
have not considered galaxies farther than 0.5 h^, 1 Mpc, to 
be consistent with the characteristics of our own data. 




0.001 - % 



Figure 7. The redshift - X-ray luminosity [0.1-2.4 keV] distri- 
bution for all clusters in the REFLEX catalogue (open circles) 
and those from the 81 cluster subsample used in this work (filled 
circles). 



5.2 Inhomogeneous redshift data sources 

The galaxy redshifts used in this work come from three 
different sources. When a cluster was already known and 
well studied, redshifts were taken from the literature, possi- 
bly from different authors too. When a cluster was already 
known but had little data available in the literature, we took 
our own redshifts in order to secure and/or improve the clus- 
ter redshift. Finally, in the case of clusters discovered in this 
project, a complete optical follow-up has been carried out. 

It is thus not unreasonable to think that such a variety 
of redshift sources might introduce some undesirable effect 
particularly on the o v estimate and therefore in the final 
L x — a v relation. To check whether this is the case, we have 
restricted ourselves to a sample of 38 clusters for which more 
than 7 redshifts per cluster were available, all of them being 
redshifts measured by us from our own data. This is then to 
be considered a completely homogeneous sample in terms of 
both luminosity and redshift, as both quantities have been 
measured in the same way for all 38 clusters. The result of 
the ODR fitting is log(L 45 ) = (-0.63 ± 0.10) + (2.2 ± 0.4) x 
log (<75oo )■ In this case the fit is quite poor, probably because 
this sample suffers from important selection effects. First, as 
we have now a small number of data points we should worry 
about the problems of measuring <j v with galaxies at small 
projected distance from the cluster centre, as explained in 
Section E3 

The minimum number of galaxy redshifts set to measure 
the a v may have a role too, specially in this case where not 
too many clusters are available for the fitting. When we 
consider the whole sample, clusters with smaller number of 
individual galaxies are given much less weight into the fit 
than those with a larger number of observed galaxies. That 
is not the case here, as all of the clusters show small numbers 



(only 3 have more than 15 measured z), so the differences 
in the weights are not so dramatic. 



5.3 Is this a representative sample from the whole 
REFLEX catalogue? 

Of course, we would like to use all 452 clusters in the RE- 
FLEX catalogue to estimate the L x — a v relation. This can- 
not be done, unfortunately, as the survey was targeted to 
measure cluster redshifts and not velocity dispersions in 
these clusters. While Multi-Object Spectroscopy observa- 
tions were used for a few distant, compact clusters (due 
to the small field of view, 5 arcmin, of the EFOSC spec- 
trograph on the ESO 3.6m telescope, see e.g. Guzzo et al. 
1999), providing around 15 redshifts per cluster, most of the 
new REFLEX cluster redshifts come from (multiple) single 
slit observations which do not deliver enough galaxy red- 
shifts to compute a velocity dispersion, although the main 
goal was always to have at least 5 spectra per target. 

We have compared the L x — z distribution for the whole 
catalogue with the one corresponding to the 171 cluster sub- 
sample used in this work. FigureQshows as filled circles the 
clusters from the subsample, and as open circles the ones 
from the whole REFLEX catalogue. One can see that both 
samples are approximately equally distributed. The results 
from a Kolmogorov-Smirnov test are compatible with this 
conclusion: we find a maximum distance value between the 
two samples (the entire catalogue and the 171 cluster sam- 
ple) of 0.11, the probability of this value being larger if com- 
ing from a random sample being 0.22 . 



The X-ray Luminosity - Velocity Dispersion Relation in the REFLEX Cluster Survey 7 



5.4 Bias from a fiux-limited sample: A 
volume-limited subsample 

A bias is expected when measuring the L x —a v relation from 
a flux-limited sample. Specifically, at the flux cut-off level 
only the brightest clusters will be considered, for a given 
velocity dispersion. As our sample has actually been drawn 
from a flux-limited sample, we may worry about whether we 
are introducing such a bias towards larger luminosities. 

This bias can be avoided by considering a volume- 
limited subsample. Therefore, we have repeated the analysis 
for the largest volume-limited sample that we could con- 
struct from the 171 cluster sample. It was built by taking 
all clusters with luminosity L x 0.5 x 10 44 h^" 2 erg s _1 , 
which corresponds to a redshift limit of z max = 0.08. To 
this redshift the comoving volume surveyed by REFLEX is 
147.4 ■ 10 6 h^, 3 Mpc 3 

We have a total of 51 "useful" clusters (according to 
the selection criteria outlined in Section |2J out of the 88 
REFLEX clusters in this volume. The ODR fitting gives 
log(L 45 ) = (-L20 ± 0.09) + (3.2 ± 0.3) x log(a 5 oo), only 
compatible at the three-sigma error level with the one ob- 
tained from the whole sample, which is an indication of some 
bias present in the total sample. 



6 THE L x - T AND a v -T RELATIONS 

Another important scaling relation in galaxy clusters is the 
X-ray luminosity-temperature L x — T. Again, measurements 
show that it does not follow the scaling relations predicted 
by pure gravitational collapse models (L x tx T 2 ). The ob- 
served L x — T relation is closer to L x oc T 3 (e.g. Novicki 
et al. 2002, Fairley et al. 2000, Xue & Wu 2000, Arnaud & 
Evrard 1999, Reichart et al. 1999), which seems to indicate 
the presence of additional non-gravitational heating sources 
in the ICM. 

o v — T has also been studied in the literature. There 
seems to be a consensus in that a v oc T~ 6 (see for exam- 
ple Xue & Wu 2000, Girardi et al. 1996, Bird, Mushotzky 
& Metzler 1995, Lubin & Bahcall 1993), again not in agree- 
ment with what would be expected in a pure gravitational 
collapse scenario, a v oc T 0,5 . 

6.1 L x — T and a v — T in the entire sample 

We have studied these relations within our cluster sample 
to see to which degree our data confirm or otherwise those 
results. ASCA temperatures from two different catalogues 
(Horner 2001 and Ikebe et al. 2002) were used for 54 clus- 
ters included in our 171 cluster sample. Ikebe et al.'s mea- 
surements were preferred whenever both works had data for 
the same cluster, as they used a two-temperature model (in- 
stead of Horner's one component model) for the isothermal 
plasma, allowing a multiphase intra-cluster medium, thus 
obtaining more accurate T estimates. 

Again, we performed the ODR fitting in the log-log 
space, obtaining the following best-fit relations (see also Fig- 
ures |H1 and GIJ : 

log(L 45 ) = (-2.53 ± 0.16) + (3.1 ± 0.2) ■ log(T) (7) 
log(a 5 oo) = (-0.45 ± 0.11) + (1.00 ± 0.16) • log(T) (8) 




kT (keV) 



Figure 8. L x — T for a subsample of 54 clusters drawn from 
the total sample, for which ASCA temperatures were available. 
The solid line is the ODR best-fit model, with one-sigma errors 
marked as dashed lines. 




kT (keV) 

Figure 9. u v — T for a subsample of 54 clusters drawn from 
the total sample, for which ASCA temperatures were available. 
The solid line is the ODR best-fit model, with one-sigma errors 
marked as dashed lines. 



The L x — T estimate that we find is in good agreement 
(at one-sigma confidence level) with the slope and intercepts 
previously measured by Fairley et al. 2000 and Arnaud & 
Evrard 1999. 

For the a v — T relation, we are finding here a larger slope 
than found by other authors (which is w 0.6). To check that 
we are not biased for including clusters with substructure in 
the sample, we select a new set of clusters from the sample 
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in Section f5. II for which temperature data are available. We 
have 33 systems, for which we find 

log(<7 5 oo) = (-0.38 ± 0.13) + (0.92 ± 0.18) ■ log(T) (9) 

in agreement with the previous result. Furthemore, we 
can select only clusters with the most reliable a v determi- 
nation, that is, those with a v measured from more than 30 
galaxy redshifts. In doing so, we end up with a sample of 16 
clusters, giving 

log((J5oo) = (-0.37 ± 0.15) + (0.9 ± 0.2) • log(T) (10) 

6.2 L x —T and a v — T in the voulme limited sample 

It is interesting to probe these relations also in a volume- 
limited sample. From the one we constructed in Sect. 15.41 we 
select those clusters for which temperature data are avail- 
able, ending up with a 28 cluster sample. On this we find 
that 

log(L 45 ) = (-2.7 ± 0.3) + (3.2 ± 0.4) ■ log(T) (11) 
log(o- 5 oo) = (0.49 ± 0.05) + (0.77 ± 0.19) ■ log(T) (12) 

The L x — T is still the same, but the a v — T is flatter than 
what is found in the 54 cluster sample, steeper than what is 
expected from self-similar models (at one-sigma confidence 
level), and in good agreement with what is found by other 
authors (Xue & Wu 2000, Girardi et al. 1996 and Lubin & 
Bahcall 1993, for example). 

In conclusion, we have derived an L x — T similar to what 
is found by previous studies, and which is in contradiction 
with self-similarity. But the a v — T we are measuring is not 
reliable, as it has shown to vary depending on the sample's 
selection criteria. The scatter and the large error bars in the 
velocity dispersions make it impossible to obtain a consistent 
relation. Still, we can draw a useful conclusion, which is that 
self-similarity is ruled out in any case. 



7 SUMMARY AND CONCLUSIONS 

We present an estimate of the L x — a v relation computed 
with a new large and homogeneous sample of clusters from 
the REFLEX survey. All studies on the L x — a v relation in 
large samples (with more than « 100 clusters) have been 
performed so far on compilations from the literature requir- 
ing the use of models to combine different data from different 
authors and/or instruments (e.g. Xue & Wu 2000, Mahdavi 
& Geller 2001). 

From the REFLEX catalogue we construct a sample of 
171 clusters for which we can derive reliable estimates of cr„, 
and then fit the log(L^) — log(er„) relation, finding 

log(L 45 ) = (-1.34±0.08)+(4.1±0.3)-Iog(ff 5 oo) ergs _1 h5 2 (13) 

in agreement with what is expected from self-similar 
models of cluster formation. 

The large size and homogeneity of the sample allows us 
to make an estimate of the intrinsic dispersion of the L x — a v 
relation. We find that Sint = 195 km s _1 . 

Measurements of L x and cr„ in clusters with multiple 
components may not be reliable. The identification in the 
literature of those clusters allows us to build a new sample 
where we remove those with significant substructure in their 



X-ray emission in Schuecker et al. (2001b). This procedure 
leaves us with a sample of 123 clusters, on which we fit again 
the log(.L x ) — \og(a v ) relation, finding 

log(L 45 ) = (-1.28±0.10) + (4.2±0.4)-log(a 50 o)ergs- 1 h5 2 (14) 

compatible at one-sigma confidence level with the slope and 
intercept that we find when using the 171 cluster sample, 
suggesting that the sample is not noticeably affected by clus- 
ters with multiple components. The fact that we are consid- 
ering only galaxies close to the cluster centre may be a reason 
for that. 

We have also investigated the L x — a v relation in a 
volume-limited sample, finding log(L 4 s) = (—1.20 ±0.09) + 
(3.2L0.3) xlog((T5oo). This significantly flatter slope is an in- 
dication that some bias may be present on the whole sample 
due to the flux limit imposed on the data. 

The L x —T and a v — T scaling relations have also been 
explored, finding that L x cx T 31±0 2 and a v cx T 100±016 , in 
a subsample of 54 clusters for which T was available from 
the literature. These results are consistent (within the error 
bars) with the one obtained for the L x — a v , as we find 
L x cx <7 3 ' 1±0 ' 7 when using these L x —T and a v — T to derive 
L x & v . 

The slope found in the L x — T relationship is steeper 
than the value predicted by a purely gravitational collapse 
model of cluster formation, and is in good agreement with 
previous measurements performed on other cluster samples. 
As for the a v — T, the slope we find is only compatible with 
the self-similar value at the three-sigma confidence level, and 
slightly larger than the value found by other authors. Also, if 
we consider a volume- limited sample, the slope found is then 
in good agreement with previous studies, yet incompatible 
with self-similarity. 

The fact that this model fails to reproduce the measured 
slope is attributed to the contribution to the ICM energy 
budget from other non-gravitational heating sources (see e.g. 
Bialek, Evrard & Mohr 2002), like AGNs, supernovae, gas 
cooling (Pearce et al. 2000), the presence of cool cores (Allen 
& Fabian 1998), or the effect of shocks in cluster formation 
(Cavaliere et al. 1997). 

Here we are faced with a paradoxical result: the slope 
of L x — a supports the self-similar scenario, whereas L x —T 
and a v — T do not, and all three relations are consistent 
with each other within one sigma confidence level. Clearly, 
more data are required to clarify the situation and reduce 
the error bars. This is a reflection of the general situation in 
this topic in the literature: so far there is a general agrement 
on L x — T and a v — T, but there is not such an agreement 
on the L x — a v . 



8 ACKNOWLEDGEMENTS 

We thank C. Altucci as her laurea thesis served as a start- 
ing point and guideline for this work. We are grateful to 
the whole REFLEX team which helped in the compila- 
tion of this catalogue. We are also indebted to Mrs. Yu 
Xing, who provided us with the data needed to compute the 
bolometric luminosity out of the luminosity in the ROSAT 
band. We also thank M. Girardi and S. Borgani for stim- 
ulating discussions and suggestions. We are indebted to T. 
Beers who gently provided us with the ROSTAT routines, 



The X-ray Luminosity - Velocity Dispersion Relation in the REFLEX Cluster Survey 9 



as well as to the authors of the ODR package, P. Boggs 
and J. Rogers. PS acknowledges support by DLR under the 
grant No. 50 OR 9708 35. This research has made use of the 
NASA/IPAC Extragalactic Database (NED) which is oper- 
ated by the Jet Propulsion Laboratory, California Institute 
of Technology, under contract with the National Aeronau- 
tics and Space Administration. This research has made use 
of NASA's Astrophysics Data System. 



REFERENCES 

Allen S.W., Fabian A., 1998, MNRAS, 297, L57 
Anders E., Grevesse N., 1989, Geochim. Cosmochim. Acta, 
53, 197 

Andrews D.F., Bickel P.J., Hampel F.R., Rogers W.H., 
Tuckey, J.W., 1972, in "Robust Estimates of Location: 
Survey and Advances" , Princeton University, Princeton 

Arnaud K.A., 1996, Astronomical Data Analysis Software 
and Systems V, eds. Jacoby G. and Barnes J., pl7, ASP 
Conf. Series, 101. 

Arnaud M., & Evrard G. E., 1999, MNRAS, 305, 631 

Beers T.C., Flynn K., Gebhardt K., 1990, AJ, 100, 32 

Bialek J. J., Evrard A. E., Mohr J. J., 2001, ApJ, 555, 597 

Bird C, Mushotzky R.F., Metzler C.A., 1995, ApJ, 453, 40 

Bohringer H. et al., 2001, A&A, 369, 826 

Bohringer H. et al., 2002, ApJ, 566, 93 

Boggs P.T., Byrd R.H., Rogers J.E., Schnabel, R.B., NIST 
IR 4834, U.S. Government Printing Office 

Borgani S., Rosati P., Tozzi P., Norman C, 1999, ApJ, 517, 
40 

Cavaliere A., Menci N., Tozzi P., 1997, ApJ, 484, L21 
Collins C.A. et al., 2000, MNRAS, 319, 939 
Danese L., de Zotti G., di Tullio G., 1980, A&A, 82, 322 
Dell' Antonio LP., Geller M.J., Fabricant D.G., 1994, AJ, 
107, 427 

Edge A.C., Stewart G.C., 1991, MNRAS, 252, 428 
Fairley B.W., Jones L.R., Scharf C, Ebeling H., Permian 

E., Horner D., Wegner, G., Malkan M., 2000, MNRAS, 

315, 669 

Girardi M., Fadda D., Giuricin G., Mardirossian F., 

Mezzetti M., Biviano A., 1996, ApJ, 457, 61 
Girardi M., Mezzetti M., 2001, ApJ, 548, 79 
Guzzo L. et al. , 1999, The Messenger, 95, 27 
Helsdon S. F.,Ponman T.J., 2000, MNRAS, 319, 933 

Horner L\ 2001 PhD thesis, 

://www-int.stsci.edu~horner/thesis.html 
Ikebe Y., Reiprich T.H., Bohringer H., Tanaka Y., Ki- 

tayama T., 2002, A&A, 383, 773 
Kaiser N. 1991, ApJ, 383, 104 
Lubin L., Bahcall N., 1993, ApJ, 415L, 17 
Mahdavi A., Bohringer H., Geller M.J., Ramella M., 2000, 

ApJ, 534, 114 
Mahdavi A., Geller M.J., 2001, ApJ, 554, L129 
Mulchaey J.S., Zabludoff A.I., 1998, ApJ, 496, 73 
Novicki M.C., Sornig M., Henry J.P, 2002, AJ, 124, 2413 
Pearce F. R., Thomas P. A., Couchman H. M. P., Edge, A. 

C, 2000, MNRAS, 317, 1029 
Peres C.B.,Fabian A. C, Edge A. C, Allen S. W., Johnstone 

R. M, White, D. A. ,1998, MNRAS 298, 416 
Quintana H ., Melnick J., 1982, AJ, 87, 972 



Reichart D.E., Castander F.J., Nichol R.C., 1999, ApJ, 
516, 1 

Schuecker P. et al., 2001a, A&A, 368, 86 

Schuecker P., B'ohringer H., Reiprich T.H., Feretti L., 

2001b, A&A, 378, 408 
Schuecker P., Bohringer H., Collins C.A., Guzzo L., 2003a, 

A&A, 398, 867 
Schuecker P., Caldwell R.R., Bohringer H., Collins C.A., 

Guzzo L., Weinberg N.N., 2003b, A&A, 402, 53 
Solinger A.B., Tucker W.H., 1972, ApJ, 175, L107 
Wainer H., Thissen D., 1976, Psychometrika, 41, 9 
White D.A., Jones C, Forman W., 1997, MNRAS, 292, 419 
Xue Y.-J., Wu X.-P, 2000, ApJ, 538, 65 



http 



